Symmetry energy and the isospin dependent equation of state 
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The isoscaling parameter a, from the fragments produced in the multifragmentation of ^*Ni + 
ssj^j^ 58pg _|_ 58j^; ^^^ 58pg _|_ 58pg j-gactions at 30, 40 and 47 MeV/nucleon, was compared with 
that predicted by the antisymmetrized molecular dynamic (AMD) calculation based on two different 
nucleon-nucleon effective forces, namely the Gogny and Gogny-AS interaction. The results show that 
the data agrees better with the choice of Gogny-AS effective interaction, resulting in a symmetry 
energy of ~18 - 20 MeV. The observed value indicates that the fragments are formed at a reduced 
density of- 0.08 fm"^. 



PACS numbers: 25.70.Pq, 25.70.Mn, 26.50.+X 



I. INTRODUCTION 

The asymmetry term in the nuclear equation of state 
(EOS) of strongly interacting matter is important for 
studying the structure, chemical compositions and evo- 
lution of neutron stars, and the dynamics of supernova 
explosion J, 2, 3J. This term dominates the pressure 
within the neutron stars at densities oi p < 2po (where 
Po is the normal matter nuclear density), and modifies 
the protoneutron star cooling rates [1, |a| . In a supernova 
explosion, the variation of the asymmetry energy can also 
change the rate of electron capture in a collapsing star, 
thereby altering the course of final explosion |y| . Under 
laboratory controlled conditions, the low density depen- 
dence of the asymmetry term in EOS can be investigated 
in a multifragmentation reaction, where the system ex- 
pands to a density of p w po/6 - Po/3 and disassembles 
into various light and heavy fragments [3,|3j!9i 10] • The 
isotopic composition of these fragments determines the 
properties of the system and has significant sensitivity to 
the density dependence of the asymmetry term. Various 
theoretical studies [13, Il3. [ij, uM have explored this in- 
fluence of the density dependence of the asymmetry term, 
however, the detailed nature of this dependence is model 
dependent. For example, the hybrid model approach pre- 
dicts that an "asy-stifP' EOS leads to fragments that are 
more neutron-rich than those produced when the EOS is 
"asy-soft" [ij]. On the other hand, calculations with the 
expanding evaporating source (EES) model predict the 
opposite trend [lH . 

Recently |lj| , the fragment yields from the nuclear col- 
lisions simulated within the antisymmetrized molecular 
dynamics (AMD) were reported to follow a scaling be- 
havior of the type. 



Y2{N,Z)/Yi{N,Z)^e'' 



N+PZ 



(1) 



where the parameters a and /3 are related to the neutron- 
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proton content of the fragmenting source. A linear rela- 
tion between the isoscaling parameter a, and the isospin 
asymmetry (Z/A)^ of the fragments, with appreciably 
different slopes, was observed for two different choices 
of the effective nucleon-nucleon interaction force. It was 
shown that such a relation can be used to limit the sym- 
metry energy parameter in the nuclear equation of state 
and extract density pertinent to the fragmentation. 

In this paper, we compare the experimentally observed 
scaling parameter a, to those predicted by two different 
effective nucleon-nucleon interaction forces, namely the 
Gogny and Gogny-AS (a modified form of the Gogny 
interaction), in the AMD calculation. It is shown that 
the Gogny-AS interaction explains the observed results 
better. The symmetry energy is found to be about 18 - 
20 MeV, indicating that the fragments are formed when 
the system expand to a density p — 0.08 fm^'^. 



II. EXPERIMENT 

The measurements were performed at the Cyclotron 
Institute of Texas A&M University (TAMU) using beams 
provided by the K500 Superconducting Cyclotron. Iso- 
topically pure beams of ^^Ni and ^*Fe at 30, 40 and 47 
MeV/nucleon were bombarded on self-supporting '"^^Ni 
(1.75 mg/cm^) and ^*Fe (2.3 mg/cm^) targets. The tar- 
gets were placed in the center of a scattering chamber 
that was housed inside the TAMU 47r neutron ball detec- 
tor jig. Six discrete particle telescopes placed inside a 
scattering chamber and centered at laboratory angles of 
10°, 44°, 72°, 100°, 128° and 148°, were used to measure 
fragments from the reactions. Each telescope consisted 
of a gas ionization chamber (IC) followed by a pair of 
silicon detectors (Si-Si) and a Csl scintillator detector, 
providing three distinct detector pairs (IC-Si, Si-Si, and 
Si-Csl) for fragment identification. The ionization cham- 
ber was of axial field design and was operated with CF4 
gas at a pressure of 50 Torr. The gaseous medium was 6 
cm thick with a typical threshold of ~ 0.5 MeV/nucleon 
for intermediate mass fragments. The silicon detectors 
had an active area of 5 cm x 5 cm and were each sub- 
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FIG. 1: Yield ratios of the fragments from Fe + Fe and 
Ni + Ni reactions (top), and Fe + Ni and Ni + Ni reac- 
tions (bottom), as a function of neutron number N for the 30 
MeV/nucleon beam energy. The solid lines are the exponen- 
tial fits to the data. The fit parameter a, is shown in the top 
left corner of each panel. 



divided into four quadrants. The first and second silicon 
detectors in the stack were 0.14 mm and 1mm thick, re- 
spectively. The dynamic energy range of the silicon pair 
was - 16 - 50 MeV for ^He and ~ 90 - 270 MeV for 
^^C. The Csl scintillator crystals that followed the sili- 
con detector pair were 2.54 cm in thickness and were read 
out by photodiodes. Good Z identification was achieved 
for fragments that punched through the IC detector and 
stopped in the first silicon detector. Fragments measured 
in the Si-Si detector pair also had good isotopic sepa- 
ration. In this work, we present data collected in the 
detector telescope placed at 44° in the laboratory. The 
fragments detected at this angle originate predominantly 
from central events. The centrality was assured by gat- 
ing on the measured neutron multiplicity. Further details 
can be found in Ref. [I^l 



III. RESULTS AND DISCUSSION 



Fig. 1 shows the yield ratios for the fragments with 
3 < Z < 7 from ^SFe + ^^Ni, ^^Pe -t- ^^Fe and ^^Ni -|- 
^^Ni reactions at 30 MeV/nucleon, as a function of neu- 
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FIG. 2: Scaling parameter a, obtained from a fit similar to 
those shown in Fig. 1, as a function of bombarding energy. 
The solid lines are exponential fit to the data. 



tron number N. The upper panel shows the ratios for the 
58pg _|_ 58pg ^j^^ 58]\jj _|_ 58]\fj reactions, and the lower 

one for the ^*Fe + ^^Ni and ^^Ni + ^^m reactions. The 
measured ratios clearly satisfy the scaling relation of the 
form shown in Eq. 1 , and are shown by the solid lines in 
the figure. The scaling parameter a is as shown in each 
panel of Fig. 1. A similar relation is also observed for 
the 40 and 47 MeV/nucleons beam energies. The values 
of the a parameters from all three beam energies for the 
two different combinations of reactions arc as shown in 
Fig. 2. The solid lines are the exponential fits to the 
data. One sees a gradual decrease in the a values as ex- 
pected due to an increase in the excitation energy or the 
temperature [iq. 

As previously discussed and shown in Ref. [i^ll3i the 
parameter a, is related to the fragment isospin asymme- 
try, (Z/A)^, through a linear relation of the form 



4C. 



sym 



T 



[{Z/A)l (Z/A) 



(2) 



where Csym is the symmetry energy and T is the tem- 
perature at which the fragments are formed. The quanti- 
ties, {Z/A)f, are the isospin asymmetries of the fragments 
in the two reaction systems, z = I and i ~ 2. Based on the 
calculations predicted in Ref. |15j , the isospin asymmetry 
of the fragments for the present systems were estimated 
at t = 300 fm/c of the dynamical evolution. The asym- 
metry values, (Z/A)^, at t = 300 fm/c for the present 
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FIG. 3: AMD calculations of the asymmetry (Z/A)^, of the 
fragments at t = 300 fm/c for the Gogny (solid line and solid 
squares) and Gogny-AS (dotted line and hollow squares) in- 
teractions at 35 MeV/nucleon. The calculations are taken 
from Ref. 1151 for the systems shown by the square symbols. 
The lines are linear fit to the square symbols. The circle sym- 
bols are the interpolated values for the systems studied in this 
work. 



FIG. 4: Scaling parameter 
of the fragments. The solid 
calculations for the Gogny 
MeV/nucleon, respectively 
interpolated values of the 
for 35 MeV/nucleon. The 
for the 30 MeV/nucleon. 
symbols. 



Q, as a function of the asymmetry 
and the dotted lines are the AMD 
and Gogny-AS interactions at 35 
|l5ll . The circle symbols are the 
parameter from the present work 
triangles are the observed values 
The error bars are of the size of 



systems were obtained by interpolating between those of 
Ref. 13 • Fig. 3 shows the AMD calculation of the frag- 
ment asymmetry, (Z/A)^ at t = 300 fm/c, as a function 
of initial asymmetry at time t = fm/c, for two different 
choices of the nucleon-nuclcon interaction force, Gogny 
and Gogny-AS. The asymmetry values for the systems 
investigated in Ref. [13 are shown by the solid and hol- 
low square symbols for the Gogny and Gogny-AS inter- 
action, respectively. The lines are the linear fits to the 
calculations. The interpolated values for the present sys- 
tems are shown by the solid and hollow circles. The solid 
and the hollow stars correspond to the reactions studied 
by Geraci et al. |20], and will be discussed later. 

It should be mentioned that the AMD calculations car- 
ried out in Ref. [l5j and shown in Fig. 3 correspond 
to the beam energy of 35 MeV/nuclcon. The interpo- 
lated values of the asymmetries for the present systems 
obtained from Fig. 3 are thus for the beam energy of 
35 MeV/nucleon. In order to make a straightforward 
comparison, we obtain the a parameter for the present 
system at 35 MeV/nucleon by interpolating the observed 
values using Fig. 2. We then explore the relation shown 
in Eq. 2 and determine the symmetry energy, knowing 
the temperature T. 

The temperatures for the present system were deter- 



mined from the double isotope ratio method of Albergo 
et al. 21] for the 30 and 40 MeV/nucleon and were ob- 
served to be 3.4 and 3.7 MeV respectively. This is con- 
sistent with the value of 3.4 MeV quoted in Ref. |l3 for 
the systems of similar mass and energy. 

Fig. 4 shows the a parameter plotted as a function of 
the difference in the asymmetry of the fragments for the 
two reactions at t = 300 fm/c. The solid and the dotted 
lines are the AMD predictions of the Gogny and Gogny- 
AS effective forces at a beam energy of 35 MeV/nucleon, 
respectively. The solid and the hollow circles are the 
interpolated values obtained from the present measure- 
ments for the 35 MeV/nucleon using two different values 
of the difference in asymmetry, resulting from the Gogny 
and Gogny-AS effective forces. For reference, the mea- 
sured values of the a parameter for the 30 MeV/nucleon 
values are also shown in the figure by the solid and hollow 
triangles. As shown in the figure, the magnitude of the 
a parameter for the 35 MeV/nucleon linearly increases 
with increasing difference in the asymmetry of the two 
systems as predicted by the AMD calculation. Also, the 
parameter values are in better agreement with the val- 
ues obtained from the Gogny-AS interaction than those 
from the usual Gogny force. A slightly lower value of 
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FIG. 5: Scaling parameter q, as a function of the asymmetry 
of the fragments. The sohd and the dotted hnes are the AMD 
calculations for the Gogny and Gogny-AS interactions at 35 
MeV/nucleon, respectively [l^. The solid and hollow circle 
symbols are the extracted values from the present work for 35 
MeV/nucleon. The various symbols are as explained in the 
text. 



the circle symbols for 35 MeV/nucleon compared to the 
Gogny-AS value (dotted hne) could be due to the sec- 
ondary de-excitation effect, which can significantly lower 
the experimentally observed value |23 • The closer agree- 
ment with the Gogny-AS type of interaction nevertheless 
indicates a stiffer EOS relative to the usual Gogny inter- 
action in the low density region of the EOS. A similar 
conclusion was also arrived at in a model dependent hy- 
brid calculations for the "^Sn + ^^^Sn, "^Sn + i24sn 
and i24sn + 1243^ reactions at 50 MeV/A Q. 

From the slope, iCsym/T (= 21.2), of the Gogny-AS 
interaction (dotted line) in Fig. 4, and temperature T ( 
= 3.4 - 3.7 MeV), we obtain a value for symmetry energy 
C^„ = 18-20 MeV. This suggests (see Fig. 1 of Ref. 
[13) that fragments are formed at a density of p ^ 0.08 
fm~'^, consistent with the fact that fragmentation occurs 
at a reduced density. A similar value of the symmetry 
energy Csym = 17-21 MeV, for the present system was 
also obtained following the statistical model calculation 
of Botvina et al. |l3 which includes the secondary de- 
excitation of the primary fragments. 



We next com par e the scaling parameter a, taken from 
various works [l9l l20l| with those obtained from the 
present study. The data were specifically chosen from 
reactions where the system temperature ranges from T 
= 3 MeV to T = 3.8 MeV, similar to those observed in 
the present study and obtained by the double isotope ra- 
tio method. The comparison is shown in Fig. 5, where 
the circle and the triangle symbols are the same as shown 
in Fig. 4 and correspond to the present study. The solid 
and the hollow stars correspond to ^^^'^^^Sn + ^'''^^Ni re- 
actions at 35 MeV/nucleon studied by Geraci et al. [23| ■ 
The solid and the hollow squares are for the H -I- ^^^'-'^^^Sn 
reactions at 6.7 GeV, whereas the inverted triangle and 
the diamond symbols are for the '^He + ^^"''^^^Sn reac- 
tions at 15.3 GeV 1^. All the data points shown corre- 
spond to reactions with nearly the same temperature as 
indicated in Fig. 5. One observes that the data points all 
lie scattered close to the prediction of Gogny-AS interac- 
tion (dotted line) . The present observation thus indicates 
that at the threshold of multifragmentation (T = 3 - 3.8 
MeV), the fragments are formed at a reduced density of 
^ 0.08 fm""^ with a symmetry energy Csym = 18-20 
MeV, and the process is consistent with an equation of 
state determined by the Gogny-AS type of interaction 
force. A more detailed study using neutron rich beams 
and populating higher asymmetry region in the a versus 
fragment asymmetry plot, where the difference between 
the predictions for the Gogny and Gogny-AS interaction 
becomes large, could be important. 



IV. CONCLUSION 

In conclusion, we have shown using experimental data 
that the linear dependence of the isoscaling parameter a, 
as a function of fragment isospin asymmetry (Z/A)^, is 
sensitive to the symmetry terms used in the nuclear equa- 
tion of state and follows a Gogny-AS interaction in the 
low density region of the EOS. Furthermore, at a tem- 
perature of T = 3 - 3.8 MeV, the fragments are formed 
at a reduced density of ~ 0.08 fm~'^ with a symmetry 
energy of the order of 18 - 20 MeV. 
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